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THE PACKING DENSITY OF THE n-DIMENSIONAL 
CROSS-POLYTOPE 

G. FEJES TOTH, F. FODORt, AND V. VfGH* 


Abstract. The packing density of the regular cross-polytope in Euclidean 
n-space is unknown except in dimensions 2 and 4 where it is 1. The only 
non-trivial upper bound is due to Gravel, Elser, and Kallus [9] who proved 
that for n = 3 the packing density of the regular octahedron is at most 1 — 
1.4... x 10“ 12 . In this paper, we prove upper bounds for the packing density 
of the n-dimensional regular cross-polytope in the case that n > 7. We use a 
modification of Blichfeldt’s method [2] due to G. Fejes Toth and W. Kuperberg 

ED- 


1. Introduction 

Let K C M” be a convex body (a compact convex set with interior points). A 
family /C = {K\, K 2 ,...} of congruent copies of K is a packing in R™ if the elements 
of /C are non-overlapping (their interiors are pairwise disjoint). The density of a 
packing is essentially the proportion of space covered by elements of the packing. 
The supremum of the densities of all packings with congruent copies of a convex 
body K is called the packing density of K and it is denoted by S(K). For a more 
detailed introduction into the basic properties of density see, for example, [8] and 

HD- 

One of the central problems of the theory of packing and covering is to determine 
the packing densities of particular convex bodies. The most important such body 
has always been the n-dimensional unit ball B n . However, the exact value of 6(B n ) 
is known only in the cases n = 2 and 3. In particular, 5(B 2 ) = 'n/y/Vl, proved by 
Time fTfiUTTl , and S(B 3 ) = n/y/l8, proved by Hales [TO]. Recently, Cohn and Kumar 
[5] proved that in dimension 24 the density of no sphere packing can exceed the 
density of the Leech lattice by more than a multiplicative factor of 1 +1.65 x 10~ 30 . 
Thus, for any practical purpose, we may consider 5{B 24 ) = 7r 12 /12L 

The current best asymptotic upper bound for S(B n ) is due to Kabatjanskii and 
Levenstem HD: 

m S(B n ) < 2( -0 - 599+o ( 1 b n _ 

For small dimensions, the bounds proved by Cohn and Elkies [4] are better but 
asymptotically their upper bound is the same as ©■ 


Date: August 8, 2016. 

fThis author was supported by the Janos Bolyai Research Scholarship of the Hungarian Acad¬ 
emy of Sciences. 

tThis research was supported by the European Union and the State of Hungary, co-financed 
by the European Social Fund in the framework of TAMOP-4.2.4.A/ 2-11/1-2012-0001 ’National 
Excellence Program’. 


1 





2 


G. FEJES TOTH, F. FODORJ, AND V. VfGH 1 


Other than the n-dimensional ball, the most interesting convex bodies are proba¬ 
bly the n-dimensional regular polytopes that exist in every dimension: the simplex, 
the cube, and the cross-polytope. The n-cube is a tile, so its packing density is 
1. However, very little is known about the packing densities of the regular simplex 
and the regular cross-polytope for n > 3. An exception is the case when n = 4, 
then the regular cross-polytope tiles R 4 (cf. Section 22 in J8]), thus its packing 
density is 1. Very recently, Gravel, Elser and Kallus [9] proved upper bounds for 
the packing density of the regular tetrahedron (1 — 2.6 ... x 10 -25 ) and the regular 
octahedron (1 — 1.4 ... x 10 -12 ). These bounds are certainly not optimal. It seems 
unclear whether the method of Gravel, Elser and Kallus can be extended to higher 
dimensions. 

We note that there has been much work done recently in order to construct 
efficient packings of regular tetrahedra, octahedra and other solids in R 3 . For 
an overview see Torquato and Jiao [19] . Dense packings of n-dimensional cross¬ 
polytopes were constructed by Rush m- Finally, we remark that, to the best of 
our knowledge, essentially nothing is known about the covering densities of the 
regular solids. 

In the next section, we will prove an upper bound for the packing density of 
the regular n-dimensional cross-polytope using known upper bounds on 5{B n ) and 
the ratio of the volumes of the cross-polytope and its insphere. Subsequently, we 
significantly improve these upper bounds for small dimensions employing a modifi¬ 
cation of the method of Blichfeldt by G. Fejes Toth and W. Kuperberg [7j. With 
these methods we establish non-trivial upper bounds on the packing density of the 
n-dimensional cross-polytope for n > 7, and we also show that the packing density 
of the cross-polytope approaches 0 exponentially fast as the dimension tends to 
infinity. 


1.1. Upper bound using the insphere. If I\ C R" is a convex body and r{K) 
is the radius of the maximum size ball contained in K, then 

Vol (K) 

< 2 > s 

This may provide a non-trivial upper bound on the packing density of a convex body 
K whose insphere is sufficiently large in volume compared to K. Using this idea, 
Torquato and Jiao [I81H9] derived upper bounds for the packing densities of the 
regular dodecahedron and the regular icosahedron and for some of the Archimedean 
solids in R 3 , see Tables III and IV in m However, it appears that they have not 
used the insphere volume ratio to investigate the packing densities of convex bodies 
in higher dimensions. 

In R 4 , the cube, the cross-polytope and the 24-cell are tiles, so their packing 
densities are all equal to 1. The insphere volume ratio method gives an upper 
bound greater than 1 for the packing density of the regular simplex. However, one 
obtains non-trivial upper bounds for the packing densities of the 120-cell and the 
600-cell using the 6(B 4 ) < 0.13126 • 7r 2 /2 bound by Cohn and Elkies [4|, see the 
numerical values in Table [l] 

In dimensions higher than 4, there exist only three regular solids, the simplex, 
the cube and the cross-polytope. The n-dimensional cube is always a tile in R ra , 
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n 

P 

Upper bound on S(P) 

4 

120-cell 

0.74972 

4 

600-cell 

0.69073 


Table 1. Upper bounds on the packing densities of 4-dimensional 
regular solids obtained from their insphere volume ratios. 


thus its packing density is 1. The insphere volume ratio method does not pro¬ 
vide a non-trivial upper bound on the density of the regular simplex in any di¬ 
mension. However, it gives an upper bound on the density of the n-dimensional 
cross-polytope, which approaches 0 exponentially fast as n tends to infinity. 

Consider the regular cross-polytope 

X n := conv (±ei,..., ±e„), 


where e*, i = 1,..., n, are the standard orthonormal basis vectors of R”. It is clear 
that 

on 

Vol (X n ) = 

n\ 

and the inradius of X n is 

r n = r(X n ) = 1 fy/n. 

We say that two non-negative sequences f(n) and g{n) are asymptotically equal 
if limn^oo f(n)/g(n) = 1. The asymptotic equality of f(n) and g(n) will be denoted 
by f(n) ~ g(n). We write f(n) <C g{n) if there exists a positive real number c such 
that f(n) < c ■ g{n) for all n. 

Consider a packing of cross-polytopes in R”. Then (jT]) , © and the Stirling 
formula yield that 


( 3 ) 

( 4 ) 


6{X n ) < Y f \^\ s{B n ) 

y ’ ~ r™Vol ( B n ) K ' 

Vol {X n ) g_0.599 ra (i +o (i)) 

— rJjVol (B n ) 

_ + ^) q — 0.59 9 ti( 1+ o(1-)) 


1 / e ' 
y/2 W2 0 198 > 
<C 0.86850”. 


Thus, we conclude that 

S(X n ) —> 0, as n —> oo 

exponentially fast. Note that the upper bound in © is asymptotic in nature and 
it says nothing about the packing density of X n in specific dimensions. In order 
to obtain concrete bounds on S(X n ), we must use specific upper bounds for 5(B n ). 
The Cohn-Elkies bounds on 5(B n ) for 5 < n < 36 (Table 3 on page 711 in [4]), 
and the almost exact value of S(B 24 ) by Cohn and Kumar [5] yield by simple 
computations the upper bounds for 6(X n ) shown in Tabled 
We will improve on these bounds in Section 3. 
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n 

Upper bound on 6(X n ) 

n 

Upper bound on S(X n ) 

24 

0.98753 

31 

0.67265 

25 

0.95416 

32 

0.63268 

26 

0.90259 

33 

0.59472 

27 

0.85275 

34 

0.55877 

28 

0.80476 

35 

0.52476 

29 

0.75871 

36 

0.49264 

30 

0.71466 




Table 2. Upper bound on S(X n ) for 24 < n < 36 using the 
insphere volume ratio. 


2. Blichfeldt’s method and its extension 


Let K C R" be a convex body. A non-negative Lebesgue measurable function 
/ : [0, oo) — > R is a Blichfeldt gauge for K if it satisfies the following conditions. 

i) In{f) ■= f R n f{\x\)dx < oo. 

ii) If {<pi : i = 1,2,...} is a set of isometries of R” such that the collec¬ 
tion {ifiiK : i = 1,2,...} is a packing, then for any x £ R” it holds that 

For technical reasons, we assume the following extra condition on /. 

iii) There exists an ?'o > 0 with f(r) = 0 for all r > r o- 

The idea of Blichfeldt [2] was that if / is a gauge for a convex body K , then 


S(K) < 


Vol (. K ) 

~UfP 


Blichfeldt applied this idea only to the unit ball in [5], He used the gauge 


/oO) 


1 — -j- for 0 < r < \/2, 
0 for r > \/2 


to show that S(B n ) < (n + 2)2 ( n + 2 )/ 2 and noted that a slight improvement of this 
bound can be obtained by the gauge 


(5) 


/» 


/ 0 (r) for r > 1, 

1 — / 0 (2 — r) for r < 1. 


For 0 < g < r(K) the inner parallel domain of K with radius g is defined as 


K- e := {x G K : gB n +xCK}. 


For x £ R™, let d(x, K_ e ) denote the Euclidean distance of x from K- Q . It is proved 
in T| that if / is a Blichfeldt gauge for B n , then for any 0 < p < r(K), 


9 e {x) = f 


d(x,K g ) 


is a Blichfeldt gauge for K. Thus, writing 
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we have 

( 6 ) 


6(K) < 


Vol (K) 

~G[bT 


This method gives an upper bound on 5(K) for each 0 < g < r(K). Our objective 
is to find, or at least estimate, the best such upper bound. 

Let K n = Vol ( B n ). For a convex body K C R" and a non-negative real number 
A, the radius A parallel domain K\ of K is the set of points in R” whose distance 
from K is at most A. The Steiner formula for the volume of K\ can be written in 
the form 

n 

Vol {K + A B n )=J2 \ n - j Kn-jV j (K), 

:=o 

where Vj(K ), j = 0,... ,n are the intrinsic volumes of K introduced by McMullen 
H3J. Note that V n (K) = Vol ( K ) is the volume of A', and 2V n -\{K) = S(K) is the 
surface volume of A. 

Let /„(/) := /(0). Using an argument that is very similar to the proof of Steiner’s 
formula, one obtains that 


(7) G(Q) = Y,e i W)Vn- j (K_ e ). 

J—o 

If 0 can be calculated or estimated explicitly, then (0 provides an upper bound 
for 8(K). We will see in the next section that it can give better results than the 
ones we can obtain from the insphere volume ratio. 

The Blichfeldt technique may be used for estimating the packing density of a 
convex body for which the intrinsic volumes of its inner parallel domain can be 
calculated explicitly or at least estimated numerically. Two classes of such bodies 
were exhibited in [7] : cylinders and the radius 1 outer parallel domain of segments. 
In the next sections, we will apply this method to the cross-polytope X n . 


3. The case of the u-dimensional cross-polytope 


Let P C R" be a polytope such that all facets of P are tangent to its insphere. 
In this case we say that P is circumscribed around its insphere. It is not difficult 
to see that if P is such a polytope, then for all 0 < g < r(P), the radius g inner 
parallel domain P~ e of P is a polytope that is similar to P with similarity ratio 
(r(P) — g)/r(P). Since the jth intrinsic volume is homogeneous of degree j, it holds 
that 

V A p - e )=(^^) 3 v An j = 0,..., n. 

Thus, the right hand side of 0 becomes a polynomial of degree n of g in the case 
that P is circumscribed around its insphere, that is, 

(8) G(g) = ^W)( n r( ^ g J V n -i(P) (0 < e < (r(P))). 

In particular, X n is circumscribed about its insphere. Betke and Henk Tj deter¬ 
mined the following formula for the jth intrinsic volume of X n : 


J + 1 


Vj + 1 

j! 




(9) 


Vj{X n ) = 2 J+1 
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where 


7 (n,j) = 


3 + 1 


o~U+ 1)* 


v^. 




dy 


i-j-i 


dx 


is the outer angle at a j-dimensional face of X n . 

Although we only defined G(g) for 0 < g < r n , it is, in fact, well-defined as a 
polynomial for all g, and thus its derivatives of all orders exist at g = 0 and g = r n . 
Elementary calculus yields that 

G\r n ) = -C- 2 In-i(m(X n ) + nr™ -1 I n (f) 

(10) = r™~ 2 (nr n I n (f) - /„-i(/)Vi(X")). 

Boroczky and Henk [3] proved that for any fixed j, 

, .n, 1 (j + 1)! (vrlnn)5 

jin, j) ~ - .. r-——, as n —> oo. 

n 1 2 VT+T n 0+l 

In particular, 

(11) Vi(V n ) ~ vWI^, as n —► oo. 

It follows from property iii) of / and the Stirling formula that 
4(/) / 0 °°/(r)r Tl ” 1 dr 


In-i(f) w n _i / 0 °° f(r)r n ~ 2 dr 

< 


k„ n f 0 °° f{r)r n 2 r 0 dr 


i-i n — 1 f Q f(r)r n ~ 2 dr 


ro 


2ne 


n 


asn-> oo, 


and thus 

( 12 ) 


4 (/) < r 0 \ —In-l(f)- 


n 


Combining (0. m and m, we obtain that 
(13) 


27re 


G'{r n ) <C r™ 2 J„_i(/) ^r 0 y —~ v^vln n 

= y/nr™- 2 I n -i(f) (r 0 V2e~ Vhin) 

< 0 


for sufficiently large n. 

Note that G(r n ) = I n (f)r™, and thus by 


6(X n ) < 


Vol ( X n ) 

Vol ( X n ) 
r"Vol (B n ) 


Vol ( B n ) 

*»(/) 


which is exactly the upper bound on 5{X n ) that we obtain using the ratio of the 
volumes of K and its insphere multiplied by the upper bound on 5(B n ) from the 
Blichfeldt gauge /. Thus, (fT31) implies that for large n the Blichfeldt method gives 
a better upper bound on 6(X n ) than the insphere volume ratio combined with the 
upper bound S(B n ) coming from the Blichfeldt gauge /. 
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4. Using the original Blichfeldt gauge function 


In this section, we will use the Blichfeldt gauge f* as defined in <[5|). Then 
/„(/*) = 1) and for n > 1, 


where 


— 


(v / 2)”(n + 1) 
Routine calculations show that 


s— 

un = — ^(V2)"(i+6 n ). 

1 -(v^-l )^ 1 (l+ ^ 


n + 1 


G'( 0) =- Ioir)V n {X n ) + / 1 (/*)K 1 _ 1 (I") 

r n 

on on i 

= —ny/n—r + 2^ 


n\ 


(n- 1)! 2 


= 0, 


and 

G"(0) = n(n ~ 1 h 0 (/*)U TO (X n ) - 2— hinVu-iiX") + 2 h{f)Vn-2{X n ) 

K T n 

Tl — 1 

= 2h{f*)Vn-2{X n ) - h(nV n -l(X n ) 

r n 

n2 n (\fn ^T ( 2\ 

= t -— 2--— arccos 1 - - 1.062097 - 1 

(n- 2)! V 2 V n) 

It is easy to check that G"(0) > 0 when n > 7. Since Vol (X ra )/G(0) = 1, our 
method provides a non-trivial upper bound on 5{X n ) in the case that n > 7. 
Furthermore, since G'{r n ) < 0, the minimum of Vol {X n ) / G{g) is attained at an 
interior point of the interval [0, r n ] and the method yields an upper bound on S(X n ) 
that is better than the one obtained from the insphere volume ratio combined with 
the Blichfeldt upper bound on 5(X n ). 

Although the quantities in d9j) cannot be calculated explicitly, they can be ap¬ 
proximated by numerical methods. By such numerical calculations, one obtains for 
G(g) a degree n polynomial in g whose maximum may be approximated (again by 
numerical methods). The results of our calculations are described in Tables [3] and 
U and they are compared to the upper bounds in Table [2] in Figure [T] 

We note that the numerical calculations suggest that the value of g at which the 
maximum of G{g) is reached tends to as n —> oo. Furthermore, it also appears 
from calculations that for sufficiently large n, the terms in which the exponent of 
g is around 2n/3 dominate the polynomial G(^=). 

Finally, we remark that we fitted and exponential function on the numerical 
results obtained from /* and got the following approximate asymptotics 

<5(X n ) < 0.87434". 


5. Concluding remarks 

The upper bounds obtained from the Blichfeldt method depend on the gauge 
function used. In Section 4, we used f* because / n (/*) can be evaluated explicitly 
for all n. However, f* is not the best possible such function, although in small 
dimensions it provides better upper bounds on S(X n ) than other known Blichfeldt 
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A 

1 - 
0.9 - 
0.8 - 
0.7 - 
0.6 - 
0.5 - 
0.4 - 
0.3 - 
0.2 - 



- 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 —► 

0 5 10 15 20 25 30 35 40 

Figure 1. Comparison of upper bounds on 6(X n ) obtained from 
insphere volume ratio (diamonds) and the Blichfeldt method with 
/* (dots) for 7 < n < 36. 


n 

Upper bound 

n 

Upper bound 

7 

0.99805 

22 

0.37264 

8 

0.98606 

23 

0.33850 

9 

0.96188 

24 

0.30697 

10 

0.92730 

25 

0.27794 

11 

0.88500 

26 

0.25129 

12 

0.83754 

27 

0.22690 

13 

0.78705 

28 

0.20462 

14 

0.73524 

29 

0.18448 

15 

0.68339 

30 

0.16586 

16 

0.63247 

31 

0.14908 

17 

0.58317 

32 

0.13398 

18 

0.53596 

33 

0.12017 

19 

0.49116 

34 

0.10770 

20 

0.44896 

35 

0.09647 

21 

0.40944 

36 

0.08635 


Table 3. Upper bounds on S(X n ) using the Blichfeldt gauge /* 
and calculated by Maplel3. 


gauges for B n . For example, it is remarked in [7] that Levenstein [12] introduced 
the following Blichfeldt gauge derived from spherical codes. Let M(n, p) denote the 
maximum number of points on S’" -1 with the property that their pairwise angular 
distances are not smaller than p. Then the following function 


fn(x) 


M (n,y>) 

0 


for 0 < 1*1 < \ / 1 2 ~ 

— I l V 1— cos<p 

for./ - —2— < 1*1. 

y 1—cosy? — 1 1 


is a Blichfeldt gauge for B n n/3 < p < tt. For a brief explanation why /„ is a 
Blichfeldt gauge, see page 726 in 7\. Kabatjanski and Levenstein (cf. Formula (52) 
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n 

density 

n 

density 

40 

5.52108 x 10~ 2 

140 

1.98099 x 10~ Y 

50 

1.72421 x 10 -2 

150 

5.36214 x 10” 8 

60 

5.19017x 10 -3 

160 

1.44520 x 10" 8 

70 

1.52250 x 10 -3 

170 

3.88033 x 10" 9 

80 

4.38143 x 10 -4 

180 

1.03837 x lO” 9 

90 

1.24242 x 10 -4 

190 

2.77031 x 10" 10 

100 

3.48295 x 10“ 5 

200 

7.37113 x 10" 11 

110 

9.66572 x 10“ 6 

250 

9.26781 x 10 -14 

120 

2.66200 x 10 -6 

500 

2.25312 x 10" 28 

130 

7.28254 x 10 -7 

750 

4.01494 x 10 -43 



1000 

6.36493 x 10“ 58 


Table 4. Upper bounds on 8{X n ) using the Blichfeldt gauge /* 
and calculated by Maplel3. 


in El) proved that 


k-\-n—2\ 

(14) M(n, ip) < _ * if cos tp < t“£, 

1 r i,fc+i 

where denotes the largest root of the Jacobi polynomial of degree k with 
parameters a = (n — 3)/2. For a definition of Jacobi polynomials see, for example, 
Formula (23) in [IT]. Using m, one can obtain Blichfeldt gauge functions for 
X n in any dimension which yield concrete upper bounds on S(X n ). However, 
these Blichfeldt gauges do not give better bounds on S(X n ) than /* up to (at 
least) dimension 300. On the other hand, in dimension 500 one obtains a better 
bound using m than with /*. We note that the calculations with da become 
computationally very demanding for higher dimensions. 

Kabatjanski and Levenstein El derived from m 

(15) M(n,ip) < (sm{<p/2))- n 2-( 0 - 599+ °( 1 » n , 

which holds for ip < 63° and is the best asymptotic upper bound on M{n,ip). We 
note that if one uses the Blichfeldt gauge f n that comes from (fl5l) . then Blichfeldt’s 
theorem yields the Kabatjanski-Levenstem upper bound ([!]) on S(B n ). Together 
with (fl3l) , this indicates that the Blichfeldt method may provide a better asymptotic 
upper bound on S(X n ) than the insphere volume ratio combined with (JT]) . 

We calculated upper bounds on S(X n ) for n < 1000 using the Blichfeldt gauge 
f n derived from (1151) (omitting the unknown o(n) term). The calculations suggest 
that the minimum of Vol ( X n )/G(g ) is attained at a value g which tends to roughly 
0.767 ... x r n as n —)• oo. We fitted an exponential function on the results based on 
which we conjecture that 

6{X n ) < 0.82886 n . 

In order to prove this we would need two things: an estimates on the o(n ) term in 
m and asymptotic formulae for the intrinsic volumes Vj(X n ) for all j = 0,..., n. 
To the best of our knowledge, no such asymptotic formulae are known at present. 
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